Domain-wall profile in the presence of anisotropic exchange interactions: 

Effective on-site anisotropy 
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Starting from a _D-dimensionaI XXZ ferromagnetic Heisenberg model in an hypercubic lattice, 
it is demonstrated that the anisotropy in the exchange coupling constant leads to a _D-dependent 
effective on-site anisotropy interaction often ignored for D > 1. As a result the effective width of 
the wall depends on the dimensionality of the system, ft is shown that the effective one-dimensional 
Hamiltonian is not the one-dimensional XXZ version as assumed in previous theoretical work. We 
derive a new expression for the wall profile that generalizes the standard Landau-Lifshitz form. Our 
results are found to be in very good agreement with earlier numerical work using the Monte Carlo 
method. Preceding theories concerning the domain wall contribution to magnetoresistance have 
considered the role of D only through the modification of the density of states in the electronic 
band structure. This Brief Report reveals that the wall profile itself contains an additional D 
dependence for the case of anisotropic exchange interactions. 



The competition between exchange and anisotropy en- 
ergies stabilizes the ground state of ferromagnetic sys- 
tems as a set of domains with different magnetizations .i 
In the region between the domains the magnetization 
smoothly changes in a way to continuously connect the 
different sectors. This configuration, known as a domain 
wall (DW), is relevant to understanding the transport 
propertiesSiSii and the response of such systems to ex- 
ternal ficlds.= However, despite the fact that much is 
known about the magnetic structure of bulk materials,— 
a framework to deal with the most general situation is 
not available;^ In some specific situations the theoreti- 
cal limitations to solve the general problem can be han- 
dled, assuming that the exchange interaction is isotropic 
while the anisotropy in the system, in general, can be 
described by an on-site, also called single-spin, interac- 
tion. With these assumptions, and using a variational ap- 
proach, Landau and LifshitzS obtained the exact form of 
a Bloch wall. They found the profile of the wall to be de- 
scribed by the expression cos9{x) = — tanh(x/A), where 
6{x) is the polar angle of a classical spin vector at posi- 
tion X and A is half the effective DW width. This kind of 
one-dimensional distribution of the magnetic moments is 
always appropriate if the surface effects can be neglected. 
Under this assumption, a system in any dimension ef- 
fectively behaves as one dimensional. This simplified 
model is in very good agreement with the experimen- 
tal results in a wide variety of investigated materials^ 
where the anisotropy mainly comes from the so-called 
crystalline field. On the other hand, the formation of 
DWs in uniaxial ferromagnets with two-spin exchange 
anisotropy has been much less studied^iS. Reference has 
mentioned in page 77 that there is no experimental ev- 
idence for the existence of this kind of anisotropy. Re- 
cently, UNiGc has shown strong indication for the pres- 
ence of this anisotropic interactionjii We can expect that 
new compounds will be discovered as part of the fast ad- 
vance that experimental magnetism is currently facing in 
the synthesis of new materials. 



The simplest Hamiltonian describing anisotropic ex- 
change interactions (exchange anisotropy) is the one- 
dimensional (ID) XXZ Heisenberg model 



^xxz — J ^'^{^i ^i-i-i + SfSfj^-^ + AS^S^^i) , (1) 



where J > and A > 1 (i.e., easy-axis anisotropy). For 
simplicity we will write Eq. ([1]) as 



ttID 

^xxz 



i+1 1 



(2) 



where Hl^^ is the one-dimensional isotropic Heisenberg 
Hamiltonian and (5 = A — 1 > 0. Recently, Yamanaka 
and Komalfl have used this model to describe the forma- 
tion of a DW in a system in any dimension. However, 
as will be demonstrated, the treatment of the anisotropic 
exchange interactions deserves a more careful attention 
and, the effective one-dimensional Hamiltonian of the 
problem for D > 1 is not the Hamiltonian ([1]). In fact, 
the effective Hamiltonian must contain an additional on- 
site anisotropy contribution whose role has been ignored 
for D > 1. This contribution has its origin in the non- 
local character of the exchange anisotropy. On physical 
grounds, it contains the information about each neighbor 
spin out of the chains that support the domain structure. 
It is our main purpose to show the correct derivation 
of the ID-effective Hamiltonian from the D-dimensional 
XXZ model and obtain a closed expression to the ef- 
fective width of the wall using both the variational^ and 
scmiclassicalS. approaches. Our results will be compared 
to those obtained in Ref. [lOl and to earlier reported 
numerical Monte Carlo ddX&M- In closing, we will con- 
sider previous theories concerning the DW contribution 
to magnetoresistance in connection to the present work. 
Our starting point will be the D-dimcnsional XXZ 
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Heisenbcrg model 



ttD _ ttD 



JS ^ S^S^ , 

(a,b> 



(3) 



where the sum goes over ah bounds between nearest 
neighbors. We will be interested in solutions with the 
boundary conditions lima^^zpoo Sa = '5'(0,0,±1), a tt 
DW in the x direction. Therefore, the magnetization 
on the sample will look the same along any spin chain 
in the y direction [two-dimensional (2D) case] or both y 
and z directions [three-dimensional (3D) case]. 

The effective one-dimensional Hamiltonian can be ob- 
tained using the fact that = {a = x, y, z) for every 
nearest neighbor of the spin in the ath site, except for the 
ones on the x axis. After summing over the sites out of 
the chains in x direction, the Hamiltonian ([3]) reduces to 
the ID version 



Hlfj = Hl^, -JSY, StSt^, J5{D -l)Y.S, 

i i 

= H]g,z-J5{D-l)Y,Sf . 



(4) 



As can be seen, the one-dimensional effective Hamilto- 
nian dH) coincides with the ID version of the XX Z model 
only in the obvious case D = 1. For D > 1, the exchange 
anisotropy leads to an additional effective on-site contri- 
bution depending on D. We notice that iJj^oj also leads 
to a similar but isotropic contribution depending on S'f 
which cannot affect the form of the wall. In the usual case 
studied in Rcfs. [s] and [o], the effective one-dimensional 
Hamiltonian is simply the ID version of the original D- 
dimensional model because the anisotropy is assumed to 
be of the on-site type. The additional on-site term in ([4]) 
makes the form and the width of the wall dependent on 
the system dimension or on the coordination number of 
the crystal structure. 

The modification of the DW profile induced by the new 
terms can be obtained following a variational approach 
as in Ref. [^. A crystalline on-site anisotropy given by 
— J5c^ Sf ^ can be added to ^ in order to consider 

i 

a more general situation. We assume small values of 5 
and 5c so that the continuum limit makes sense. Then, 
Eq. HI with the crystalline on-site term added, defines the 
one-dimensional variational problem given by 



T re 

-{d^Sf + —{d^S'f - J5t SI '^]dx^ min , (5) 



where St = 5 D + Sc is the total anisotropy parameter 
depending on D and x is in units of the lattice constant. 

As it is shown in ([5]), the exchange is anisotropic and 
the on-site interaction depends on the system dimension. 
Therefore, the DW profile and its effective width are dif- 
ferent from those obtained in Ref. [s]. As mentioned be- 
fore, the on-site term reflects the fact that each site effec- 
tively feels the cost of energy to have its nearest neighbors 



pointing in directions away from the anisotropy axis. As 
the number of nearest neighbors increases with the di- 
mension we can expect the effective width of the wall to 
decrease for higher D. 

Using the angle representation for the magnetic mo- 
ments 



S{x) = S [0,sme{x),coa9{x)] , (6) 
where we assume 3> 1, Eq. ([5]) becomes 



1 



'l + (5sin2 61) 0'^ -5t cos^i 



dx = min , (7) 



which generates the equation of motion 



2(5. 



* sin 2e+{l + S sin^ 9) 0" = 0. (8) 



Integration of ([8]) with the boundary conditions of a tt- 
DW leads to the solution for the magnetic structure of 
the form 



- \/2{S D + Sc) X = arctanh 



cosw X 



^Jl + S sin^ 0{x) 



arctan 



\/S cos0{x) 
_^J^+S sin^0ix) 



(9) 



As can be seen, the standard Landau and LifshitzS ex- 
pression cos 0{x) = — tanh{\/26c x) can be obtained from 
Eq. ([9]) for (5 = 0. Equation ([9]) is thus the generalization 
of that expression for the case of anisotropic exchange 
interactions. 

The same expression for the DW profile can be ob- 
tained following a different approach, as in Ref. d. In 
that case, the static solution of the spin quantum equa- 
tions of motion, 



dSl 
dt 







(10) 



is considered. Here H can be either the ZJ-dimensional 
Hamiltonian ([3]) or the effective one-dimensional Hamil- 
tonian Qj. The sum over the nearest neighbors of the 
ath spin resulting from (llOft can be transformed in the 
continuum limit as 



Y,S^^dlS'^ix) + 2DS'^{x). 



(11) 



For the case of the effective one-dimensional Hamiltonian 
([4]), D must be taken equal to 1 in Eq. (fTTI) . Then 



the classical vector representation ([6]) can be used and 
Eq. ([H]) -and thus Eq. 1^- is recovered. We note that 
according to Eq. ([6]) one finally takes S''^{x) = 0. Thus, 
only the evolution of actually needs to be considered 
in Eq. (fTO|) . 

We can extract the effective width of the wall {W) from 
the behavior of Eq. ([9]) near the origin. For a; — > we 
get + 6 cos6{x) = —yj2{5D + 5c) x and, therefore, 



100 



W = 2 



2St 



2{6D + Sc) 



(12) 



As expected, an increase in the dimension results in a 
narrower DW. This is in agreement with the fact that 
the local field produced by the spins out of a given chain 
increases the energy cost of the spin twist. As a conse- 
quence, the system effectively gets less spins out of the 
z axis and thus the effective width of the wall becomes 
smaller. The specific behavior cx 1/y/St agrees with the 
limit case 5 — > 0, in which the anisotropic exchange term 
{J5/2){dxS^y' in Eq. ^ can be neglected. In that case, 
one gets an isotropic exchange parameter J/2 and an on- 
site anisotropy with the effective parameter J St- The 
ratio of these two parameters determines^ the effective 
width of the wall. Slightly away from that limit, Eq. 
(I12p shows an additional exchange-type dependence on 
5 (i.e., 5 is in the numerator) coming from the total ex- 
change parameter 1 + 5 = A of the z components in the 
XXZ model. 

When (5 = 0, Eq. ([12]) gives 2 ^^^ , which has no de- 
pendence on D and is the effective width of the wall in the 
well-known case for which the anisotropy in the original 
model is of the on-site type.=>2 As can be seen from fl!2|) . 
the effective width has less sensitivity to variations in the 
parameter 5 than in Sc- This is a consequence of the com- 
bined role of 6 as exchange and anisotropy parameter at 
the same time, as can be seen from the original model 
®. 

Now we are going to compare our results for D = 1 
with those obtained in Ref. [KJ, where (5c = was consid- 
ered. Following a pure quantum method in treating the 
XXZ model ([T]), Yamanaka and Komatf rederived the 
Landau and Lifshitz expression cos9{x) = — tanh(x/A) 
with an effective wall width given by 



2A = -2- 



1 

In q ' 



(13) 



where < g < 1 and q+l/q~ 2 A. 

As discussed above, we can get such an expression for 
the domain wall profile, which corresponds to the case 
of isotropic exchange interactions and on-site anisotropy, 
only in the limit case (5^0. Indeed, it is not hard to see 
that for D = 1 and i5 ^ 0, Eq. ([9]) -with 5c = 0- reduces 
to that expression with A given by 1/ ^/25. 

At first sight the DW width (fT3D looks quite different 
from our result given by Eq. (fT2D . with 5c = Q-, in the 
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Fig. 1: The effective DW width, as obtained from Eq. ([14]) 
for D = Z (open triangles), is compared to Monte Carlo data 
(filled squares) from Ref. [l2|. Open squares correspond to the 
high anisotropy limit for which Eq. (|14p is not applicable. 
Please see text for further details. 



_D = 1 case. However, as can be seen they also coincide 
for small values of (5 (Aw 1). One first notices that since 
A > 1, then (7 = A — VA^ — 1. We can expand Eqs. 
([12]) and ([13]) in powers of 5 and obtain for both cases 
2/^/25^0{5^/'^)M.T\ms, up to leading order in 1/(5, both 
the expression ([12]). with 5c = 0, and the one reported 
in Ref. [l^ are equal to 2-^^. It is interesting to note 

that the Eq. ([9]) for the form of the wall and ([12]) for the 
effective DW width have been derived for large values of 
spin S. However, both of them reduce, for Z3 = 1 and 
small values of the parameter 5, to the results of Ref. [l^ 
for spin 5 = 1/2. 

The anisotropic Heisenberg model has been studied be- 
fore using the Monte Carlo method.= Recently, Serena 
and Costa-Kramer= have made use of this technique to 
describe a tt DW in the XXZ model. In order to com- 
pare our results with theirs, we first make the equivalence 
between our parameters and those used in Ref. [l2]. Here- 
after, to avoid confusion with A in the present work, we 
identify the anisotropy parameter of Ref. [l2| by A'. One 
can see that they are related by A = 1/(1 — A'). There- 
fore (5 = A - 1 = A'/(l - A') and, from Eq. ([12]), with 
i5c = 0, we get 



W 



1 



^J2M D 



(14) 



Here, the factor 2 of Eq. ([12]) has been changed to tt 
to follow the definition used in Ref. [lH for the effective 
DW width. Equation p4p gives IF in a simpler manner 
than the one given by Eq. ([12]) for 5c = 0. Addition- 
ally, it looks closer to the result '''-^55= of the usual case 
of isotropic exchange interactions and on-site anisotropy. 
Equation ([T4]) confirms the; power-law W cx A'~^/^ that 
has been suggested in Ref. [l2|. It also elucidates the fact 
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that the proportionahty factor of this behavior is deter- 
mined by the dimensionahty D of the system. 

Figure [U compares the effective DW width, as obtained 
from Eq. (fT4|) . to the numerical data from Ref. [l2| for 
the case D ~ 3. Here some remarks arc in order. Equa- 
tion p4p has been derived in the contimmm hmit, i.e., 
large values of W and thus small anisotropy values have 
been assumed. We need A' <C 1/(21?) in Eq. ([H]) to 
get W ^ 1. Therefore, for Z? = 3, we can expect Eq. 

to give good results for A' ^ 1/6 « 0.2. This agrees 
with the fact observed by Serena and Costa-Kramer that 
for anisotropy values larger than A' = 0.5, the magne- 
tization suddenly changes in one lattice constant .= The 
value W = TT has been defined in Ref. [l2| for this high 
anisotropy limit. Open squares in Fig. [1] correspond to 
this case and are not considered in Eq. (I14p . An unex- 
pected departure from Eq. ([14]) is observed for Monte 
Carlo (MC) data with very small anisotropy values. Ac- 
cording to Ref. [m, this may be associated to numeri- 
cal uncertainties for such very low values of anisotropy 
(A' < 0.001). With these comments in mind, we con- 
clude that our results are in very good agreement with 
the numerical data of Serena and Costa-Kramer. 

Previous theories have considered the role of a DW in 
the transport properties of ferromagnetic metals. The ef- 
fective width W, when compared to the Fermi wavelength 
Af, determines how strong the effect of the DW is in 
increasing^!^ or decreasing^ the electrical resistance. Re- 
cently, Sil and EntellS have considered the modification 
of the electronic band structure by the DW and shown 
that there can be an increase or decrease of the resistance 
depending on the dimensionality, position of the Fermi 
energy and strength of the coupling between conduction 
electrons and DW. The role of D has been considered 
through the modification of the density of states. One 
can see that these theories have considered no depen- 
dence of the effective DW width on the dimensionality 
D of the system. Here we have seen that this is valid 



only for the usual case of isotropic exchange interactions 
and on-site anisotropy. This Brief Report reveals that 
the wall profile itself contains an additional D depen- 
dence for the case of anisotropic exchange interactions. 
In this case Eq. ([12]) will define such a D dependence of 
W, which is the relevant DW parameter in the electronic 
transport. 

To summarize, we have followed both the variational 
and semiclassical approaches to describe a tt DW in the 
XXZ model in a ZJ-dimensional system. It was shown 
that the effective one-dimensional Hamiltonian contains 
an on-site anisotropy contribution which has been previ- 
ously ignored for D > I. The standard Landau-Lifshitz 
DW profile for isotropic exchange interactions has been 
generalized. Increasing the dimension of the system de- 
creases the width of the wall. This case is thus essen- 
tially different from the well known one for which the 
anisotropy in the original model is only of the on-site 
type. Very good agreement with earlier numerical data 
was found. This work adds a new source of D dependence 
for measurable physical properties like the magnetoresis- 
tance. 
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